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in 

The fractional moment method, which was initially developed in the discrete con- 
text for the analysis of the localization properties of lattice random operators, is 
^ 1 extended to apply to random Schrodinger operators in the continuum. One of 

QQ ' the new results for continuum operators are exponentially decaying bounds for 

the mean value of transition amplitudes, for energies throughout the localization 
C ~ 1 , regime. An obstacle which up to now prevented an extension of this method to 

• the continuum is the lack of a uniform bound on the Lifshitz-Krein spectral shift 

f^*^ ' associated with the local potential terms. This difficulty is resolved through an 

analysis of the resonance-diffusing effects of the disorder. 

o ' 

1. Introduction 

i-^H ' The addition of disorder through a random potential may have a drastic effect on the spectral 

and dynamical properties of a Schrodinger operator. In certain energy regimes the spectrum 
of the operator may turn from absolutely continuous into dense pure point with localized 
eigenstates. This phenomenon, known as Anderson localization, also manifests itself in the 
form of dynamical localization, that is the non-spreading of wave packets supported in the 
corresponding energy regimes. 

There are two known approaches to the mathematical analysis of localization properties 
for multidimensional random Schrodinger operators. Both were initially developed in the 
discrete context, i.e. for random lattice operators. The method of multiscale analysis goes 
back to the ground breaking work of Frohlich and Spencer [5] from 1983 and, by now, has 
lead to a multitude of results on spectral and dynamical localization for a wide range of 
models. In 1993, Aizenman and Molchanov [2] introduced the fractional moment method'mto 
the study of Anderson localization. For discrete systems this method has provided a simple 
perspective on localization and has enabled exponentially decaying bounds on expectation 
values of various propagation kernels. 

Multiscale analysis has meanwhile been extended to continuum Anderson- type models, 
see e.g. [7, 4] and, for a state of the art account, [6]. For an introduction to multiscale 
analysis (which is not used in our work) and many references, see also the recent book [10]. 

Our goal here is to outline a continuum version of the fractional moment method and 
its consequences, with results roughly corresponding to those obtained for the lattice case 
in [3]. We focus on a continuum Anderson- type model in L 2 (M. d ) of the form 

H U :=H + \V U . (1.1) 
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In Section 2 we will introduce a prototypical set of assumptions on the background 
operator Ho and random potential Vu,. We then state three results: First, in Section 3, we 
discuss a crucial boundedness result for fractional moments of "smeared" Green functions, 

1. e. operator norms of spatially localized resolvents. Exponentially decaying bounds on this 
quantity will then be shown to imply pure point spectrum with exponentially decaying 
eigenfunctions as well as exponential decay for the mean values of transition amplitudes 
(Section 4). A finite volume criterion for exponential Green function bounds is provided in 
Section 5. Finally, as sketched in Section 6, applications of the fractional moment method 
are found by verifying the finite volume criterion in the usual large disorder, Lifshitz tail or 
band edge regimes. 

A full account of this work with detailed proofs of all the results stated below is provided 
in [1]. For a discussion of other consequences of the localization results established here, 
e.g. Kubo conductance and quantum Hall effect, see the contribution of A. Elgart to this 
volume. 

2. A prototypical model 

The operator H may incorporate deterministic magnetic and electric potentials, i.e. have 
the form 

H Q = ( i \7-A(q)) 2 + V Q (q). (2.1) 

Simple and for our considerations suitable assumptions are local boundedness of the vector 
potential A, its derivatives diA and the positive part Vb.+ of the electric potential. We 
also assume that Vo is bounded from below. Thus H is bounded below and we let E := 
inf a (Ho). 

Disorder is introduced into (1.1) through a parameter A and an Anderson- type random 
potential 

K,(«) = $>a iU tfa(«) ■ (2-2) 

For simplicity we will assume here that X = Z d and that U a (q) =U(q — a), a £ X, for a 
non-negative, bounded and compactly supported single site potential U, say supp{7 C Bq, 
where B r x = {q : \q — x\ < r}. For technical reasons we also assume that |d(supp£/)| = 
(dA denoting the boundary of a set A and | • | d-dimcnsional Lebesgue measure) and that 
the U a cover space in the sense that 

< b- < U c(l) < b + < 00 (2-3) 

uniformly in q £ M. d . 

Finally, we assume that the coupling parameters r) a , a £ X, are independent, identically 
distributed random variables with absolutely continuous distribution dn(rj) = p(r])dr). The 
density p is bounded and supported in [0,1]. 

A number of these assumptions can be weakened, in particular those on X, U a and r\ a . 
The coefficients of the background operator Ho may include certain L p -type singularities. 
For more discussion on what is technically necessary sec [1]. 
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3. Finiteness of fractional moments 

A central object in the fractional moment approach to localization for lattice operators is 
given by the fractional moments E(\GE+io( x ,y)\ s ), where G z (x,y) = (5 X \(H — z) _1 |<5j,) is 
the Green function, < s < 1, and E denotes averaging over the disorder. Finiteness of 
the fractional moments is seen relatively easily for suitable distribution of the random pa- 
rameters as the singularities of the Green function become integrable through the exponent 
s < 1. 

As noted previously in the multiscale analysis approach, for continuum models a useful 
counterpart of the discrete Green function \Ge+io(x, y)\ need not be the integral kernel of 
(H — E — iO) -1 , but rather the operator norm ||xa;(i? — E — iO) _1 Xj/|| for suitable compactly 
supported functions \x and \y Finiteness of their fractional moments is a crucial prelim- 
inary result for our discussion and technically deeper than the corresponding result in the 
discrete case. 

To state this result, for an open set SI C ff 1 we denote by the restriction of H = H w 
to L 2 (Q) with Dirichlet boundary conditions. Throughout we also denote characteristic 
functions of a set A by 1a and, for x E K d , \x = ^b^, where the size r of the bumps U a 
serves as a convenient length scale. 

Lemma 1. Let H u be a random Schrodinger operator as in (1.1) with assumptions as in 
Section 2. Then, for each s E (0, 1) and A > 7 there exists C s ^\ < oo such that 

supE (\\Xx H(n) ^ E _ - £ XX) < C S , A (1 + l/\Y(l + \E- £ |) s(d+2) (3.1) 

for any open f2 C R d , x, y € fi and E E R. One can choose 

C s A < const } \ ' . (3.2) 

1 — s 

As in the discrete case, the proof of this result proceeds by showing that the independent 
variation of some of the random parameters n a resolves singularities which are due to the 
proximity of the given energy to an eigenvalue whose eigenvector has significant support 
nearby. However, a change in a parameter can also have the opposite effect, through the 
creation of a resonance. In the discrete setup the latter possibility occurs at not more than 
a single value of the random parameter, since each coefficient affects a rank-one term and 
the number of energy levels which can be moved past E is bounded by the rank of the 
perturbation. Aside from the fact that the rank-one analysis is not applicable, the source 
of the difficulty in the extension of the previous analysis can be traced to the fact that, in 
the continuum setup, there is no uniform bound on the corresponding "spectral shift" . To 
circumvent these difficulties we employ the Birman-Schwinger principle in place of rank- 
one analysis, and control the Lebesgue measure of the nearly-singular values of a coupling 
parameter by means of the following "weak 1-1" type bound 

|{»j : \\T{n + A + *0)- 1 T|| HS > t}\ < j\\T\\ 2 HS , (3.3) 

valid for any maximally dissipative operator A and Hilbert-Schmidt operator T. This result 
was proven in [8]. 
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As a consequence of this analysis we find that it suffices to average over "local environ- 
ments" of x and y. Rather than taking the full expectation one merely averages over the r\ a 
with a in suitable neighborhoods of x and y. This yields a bound as in (3.1) with constants 
which arc uniform in the values of the remaining random parameters, an improvement of 
Lemma 1 which is important in the proof of Theorem 2 below. 

4. Localization properties 

The uniform fractional moment bound (3.1) holds for all energies in the class of continuum 
Anderson models considered here. In the following we will identify the existence of expo- 
nentially decaying bounds (in \x — y\) for the left hand side of (3.1) as a characteristic of 
the localization regime. We first show that such bounds for finite volume operators (but 
uniformly in the volume) imply spectral and dynamical localization. 

Let C M. d be open and A„ C fi, n G N, a sequence of bounded open domains such that 
U A„ = Q and ij( A ™) converges to ij( n ) in strong resolvent sense. We also define a modified 
distance by 

dist n{x, y) := min{|a; - y\, dist (x, tt c ) + dist (y, Q c )} . (4.1) 
Pj(H) denotes the spectral projection onto J for H and j| • ||^ r the trace norm. 

Theorem 1. Let H , and A n , n G N, be as above. Suppose that for some < s < 1 and 
an open bounded interval J there are constants A < oo and ji > such that 

j/ [}^ H ^)-E Xv¥ ) dE ~ Ae-» di8t ^*rt (4.2) 
for all n 6 N and x,y G A„. Then for every r < 1/(2 — s) there exists A r < oo such that 

E[ sup \\ Xx g(HM)Pj(HW)xy\\tr) < A r e~^ dist ^ (4.3) 
\g--\g\<i J 

for every x,y G O. Here the supremum is taken over all Borel measurable functions g which 
satisfy \g\ < 1 pointwise. 

In the case ft — R d it further holds that the spectrum of H in J is almost surely pure 
point, with eigenf unctions tp which for every v G (0, 2/(2 — s)) satisfy 

limsup h 1 MM<_,. (44) 

\x\^oc \ x \ 

The bound (4.3) with g(H ) = e ltH implies dynamical localization with exponential decay 
of the transition amplitudes for wave packets with energies restricted to J . This is stronger 
than the dynamical bounds which can be obtained through the multiscale analysis approach, 
e.g. [6] for the best known result. 

Theorem 1 as well as Theorem 2 below are applicable even when the operator exhibits 
extended boundary states in certain geometries, provided there is "localization in the bulk" . 
This is the relevance of the domain adapted metric distn- Note that distRd(a;,y) = \x — y\. 

While typical applications of Theorem 1 (see Section 5) will work with exponential 
bounds for E(\\x x (H^ - E^XyW 3 ) which are uniform in E G J , it is interesting to note 
that Theorem 1 only requires the energy-averaged bound (4.2). 
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The proof of Theorem 1 in [1] proceeds by first verifying the bound (4.3) for the fi- 
nite volume operators H^- An \ with constants uniform in n. In finite volume the norm of 
Xx9{H)Pj(H)x y may be estimated in terms of sums of bounds on rank-one operators. The 
latter have equal operator and trace norms, which ultimately allows to state (4.3) as a trace 
norm bound. 



5. A finite volume criterion 

In applications of Theorem 1 it is necessary to find energy regimes in which the exponential 
resolvent bound (4.2) can be verified. In this section we provide a finite volume sufficiency 
criterion for the desired exponential decay. 
We define the boundary layer of a set 

SA:={q:r< dist (q, A c ) < 23r} , (5.1) 

where the choice of the depth is somewhat arbitrary, but convenient for the technical im- 
plementation of the proof of the following result. 

Theorem 2. Let H be as above. Then for each s E (0, 1/3) and A > there exists M s ^\ < 
co, such that if for some E eK and L > 24r ; 

e-'< := M S;A (1 + 1/A) 2s (l + \E - E \f s(d+2) {1 + L) 2 ^ 



x limsupsupE \\ Xa TTlRL , ~Ub^ 

< 1 , (5.2) 
then for any open ft C M. d and any x,y € Q 

limsupE (||Xx g(a) _ 1 E _ . e X„ir) < e^A(s,X,E)e^ dist ^^ 2L , (5.3) 

with A(s, A, E) the right hand side of (3.1). One may choose 

(1 + A) 5s ( d + 4 ) 

M S , A = const 1 \ ' . (5.4) 

1 — 3s 

The proof of Theorem 2 proceeds essentially by an iterative argument where the distance 
from x to y is covered by balls of radius L. The bound (5.2) serves as an initial decay 
estimate for the resolvent between the center and boundary of a ball of radius L, reflected 
in the positive exponent 7. An iterative geometric resolvent expansion is used to show that 
the decay adds up (or better: multiplies up) to exponential decay with rate proportional to 
7 and dist n {x, y). As this resolvent expansion does not work near the boundary of ft, one 
uses the modified disto. 

The factors which appear in each step of the resolvent expansion are not all independent. 
Therefore a (triple) Holder bound is used to factorize their expectations. This is the reason 
for having to work with s < 1/3. In order to not having to divide the exponent s by three in 
each step of the iteration (which would cause it to collapse into 0), the random parameters 
7] a near the boundaries of domains used in the expansion are re-sampled in each step of the 
iteration. This means that they are replaced with parameters rj a which are independent 
of the 7] a , but have the same distribution, a procedure which also appears, for example, in 
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the spectral averaging argument of [9]. This allows to avoid the use of various versions of 
"decoupling lemmas" which have entered the fractional moment method for lattice models 
and seem to be harder to verify in the continuum. 

As opposed to the use of an iteratively increasing sequence of length scales in the mul- 
tiscale analysis approach, only one length scale L is used by the fractional moment method 
to go from finite to infinite volume. In this iteration process, Lemma 1 plays a role similar 
to Wegner estimates in multiscale analysis. It provides a worst case bound on the growth 
of the resolvent over distances less than L, where (5.2) can not yet be used. 

The exponential decay bound (5.3) on resolvents is not only a necessary consequence of 
the finite volume criterion (5.2), but for ft = R d is also sufficient for it, as shown by the 
following result: 

Theorem 3. Let H be as above and suppose that for some A < oo, /i > and £el 



for all a, (3 G Z d . Then, for sufficiently large L, (5.2) is satisfied uniformly for all E' in an 
open neighborhood of E. 

A particular consequence of this is that Theorem 1 could be stated under assuming the 
infinite volume exponential decay bound (5.5) for the resolvent, as (5.5) implies (5.2) and 
(5.2) implies (5.3) (which allows for finite volume) and thus (4.2) for a neighborhood J of 
E and all A„. 

It is interesting to note that Theorem 3 allows to conclude localization on an open 
interval from a bound for a single energy. This is due to the fact that for finite volume A 
the fractional moments E(\\x x (H^ — z)~ 1 x y \\ s ) arc Holder continuous in z G C. Thus the 
set of energies where a bound like (5.2) is valid must be open. 

6. Applications 

Applications of our method consist in verifying the bound (5.2) in concrete energy regimes. 
In this sense (5.2) is a fractional moment version of the initial length bounds used to start 
a multiscale analysis. Here are examples of regimes where (5.2) can be verified (for detailed 
statements see [1]): 

• The band edge/Lifshitz tail regime: Here (5.2) follows from smallness of the density 
of states in a suitable energy interval. One may work with smallness of the expected 
number of eigenvalues of finite volume operators H ( A ) or directly with smallness of 
the integrated density of states in infinite volume, e.g. Lifshitz tails. 

• The large disorder regime: Under somewhat stronger assumptions on the distribu- 
tion of the T] a , for example in the case of uniform distribution on [0, 1], one can 
improve the bounds (3.2) and (5.4) and obtain C Sj a and M Sj > which are bounded as 
A — > oo. This in turn may be used to prove localization in the large disorder regime: 
For every E' G M there exists A' sufficiently large such that for A > A' the energy 
interval (—oo,E') is localized (i.e. (5.2) can be verified for all E G (— oo,E')). 

• The multiscale analysis regime: One may also use the typical output of a multiscale 
analysis to verify (5.2) at sufficiently large L. Thus one gets the stronger dynamical 




(5.5) 
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localization bounds provided by the fractional moment method in all regimes where 
a multiscale analysis can be carried through and our general setup from Section 2 
holds. 
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